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ABSTRACT 

We study the evolution of disk accretion during the merger of supermassive 
black hole binaries in galactic nuclei. In hierarchical galaxy formation models, 
the most common binaries are likely to arise from minor galactic mergers, and 
have unequal mass black holes. Once such a binary becomes embedded in an 
accretion disk at a separation a ~ 0.1 pc, the merger proceeds in two distinct 
phases. During the first phase, the loss of orbital angular momentum to the 
gaseous disk shrinks the binary on a timescale of ~ 10 7 yr. The accretion rate 
onto the primary black hole is not increased, and can be substantially reduced, 
during this disk-driven migration. At smaller separations, gravitational radiation 
becomes the dominant angular momentum loss process, and any gas trapped 
inside the orbit of the secondary is driven inwards by the inspiralling black hole. 
The implied accretion rate just prior to coalescence exceeds the Eddington limit, 
so the final merger is likely to occur within a common envelope formed from 
the disrupted inner disk, and be accompanied by high velocity (~ 10 4 kms -1 ) 
outflows. 

Subject headings: accretion, accretion disks — black hole physics — gravitational 
waves — galaxies: active — galaxies: nuclei — quasars: general 
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1. Introduction 

A large fraction of galaxies in the local Universe harbor supermassive black holes (Magor- 
rian et al. 1998), with the black hole mass correlating well with the velocity dispersion of the 
galaxy (Ferrarese & Merritt 2000; Gebhardt et al. 2000). One explanation for this correlation 
is that black holes, like galaxies themselves, have grown in mass through the hierarchical 
merging of smaller progenitors (Haehnelt & Kauffmann 2000; Cavaliere & Vittorini 2000; 
Menou, Haiman & Narayanan 2001). This implies that the frequent galactic mergers at high 
redshift could also involve the inspiral, and probably eventual coalescence, of supermassive 
black holes. Typically, these binaries will pair black holes of disparate masses, since minor 
galactic mergers are more common than major mergers in hierarchical models of galaxy 
formation (Haehnelt & Kauffmann 2000). 

For a black hole binary with masses of 10 8 M & and 10 6 M & , gravitational radiation will 
lead to merger within a Hubble time if the separation is < 10~ 2 pc. Reaching such a small 
separation, following a galactic merger, requires that the black holes lose almost all of their 
orbital angular momentum to stars, dark matter, or gas. Angular momentum loss to stars 
occurs readily at large radii, but slows down at pc-scale separations if the supply of stars 
on low angular momentum orbits, which can interact with the binary, becomes depleted 
(Begelman, Blandford & Rees 1980; Mikkola & Valtonen 1992; Makino 1997; Quinlan & 
Hernquist 1997). The extent of this loss cone depletion remains a subject of active debate 
(Milosavljevic & Merritt 2001; Yu 2002; Zhao, Haehnelt & Rees 2002), with recent work 
suggesting that stellar dynamical processes may suffice to bring binaries into the gravitational 
radiation dominated regime. Irrespective of these uncertainties, gas, if it is present at small 
radii in the nucleus, is likely to be the most efficient sink of binary angular momentum 
(Begelman, Blandford & Rees 1980; Pringle 1991; Ivanov, Papaloizou & Polnarev 1999; 
Gould & Rix 2000). In this letter, we study the role of gas in hastening the merger process, 
and the influence of the binary on the accretion rate and luminosity of the black holes. 



2. Timescales for black hole mergers 

When the mass ratio of the binary, q = M 2 /Mi <C 1, the interaction of the lower mass 
hole with the disk is identical to that of a Jupiter mass planet in a protoplanetary disk 
(ignoring relativistic effects due to black hole spin). This regime, which Gould & Rix (2000) 
described as 'planet-like', has been extensively studied (Goldreich & Tremaine 1980; Lin & 
Papaloizou 1986). The secondary black hole will open a gap in the disk around the primary, 
and exchange angular momentum with the gas via gravitational torques. If the disk mass 
is at least of the same order of magnitude as the secondary, orbital migration then occurs 
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on the viscous timescale of the accretion disk (Syer & Clarke 1995; Ivanov, Papaloizou & 
Polnarev 1999). The radial velocity of fluid in the disk is given by v r = — 3u/2r, where v is 
the 'viscosity' used to approximate angular momentum transport within the disk at radius 
r. Using the Shakura & Sunyaev (1973) a prescription, v = ac s h, where c s is the disk sound 
speed and h the vertical scale height, the rate of migration is approximately 
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where vk is the Keplerian velocity. Equating this to the rate of inspiral due to gravitational 
radiation, which for circular orbits is 

64G 3 M 1 M 2 (M 1 + M 2 ) 
° GW = 5^ ' (2) 

we obtain the critical semi-major axis within which gravitational radiation is the dominant 
process leading to orbital decay. For g< 1 we find 
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This transition radius between disk-driven migration and gravitational wave inspiral scales 
with the gravitational radius of the primary black hole, and decreases for more extreme 
mass ratios. For fiducial parameters of (h/r) ~ 10~ 2 , a = 10 -2 , and q = 10~ 2 , we obtain 
Ocrit ~ 10 2 GMi/c 2 , or around 3.5 x 10~ 3 pc for a primary black hole mass of 5 x 10 8 M . 

Once gravitational radiation becomes significant, the dew ^ a ~ 3 dependence of equation 
(2) ensures that disk torques rapidly become negligible. The merger timescale from a cr i t is 

w^(fHr /5 «M^)> 

which is around 10 5 yr for the parameters above. During this final stage of orbital decay, 
gravitational torques will still extract angular momentum from the disk interior to the orbit, 
maintaining a gap and forcing the inner disk to accrete more rapidly. Simulations show 
that gap opening occurs on an orbital timescale, so there is no dynamical reason why this 
driven inflow should not continue until the inspiral velocity reaches a significant fraction of 
the Keplerian velocity. Adopting the criteria that dew ^ W~ 2 vk, a gap could remain until 
a ~ 17g 2//5 (GMi/c 2 ) <^ a cr i t . In practice, the thin disk approximation will instead break 
down when the accretion rate through the inner disk exceeds the Eddington limit. 

For mass ratios q > 1CT 2 , a crit is substantially outside the marginally stable orbit. 
However, we note that for smaller q, disk-driven migration may always be more important 
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than gravitational radiation. This means that, if a geometrically thin accretion disk is 
present, low frequency gravitational waves would not be detectable from extreme mass ratio 
binaries, except during the final dynamical stage of the merger. Geometrically thick accretion 
flows, however, are expected to be much less efficient at driving migration (Narayan 2000). 



3. Numerical merger model 

We model the migration of the secondary black hole in the planet-like regime (Gould 
& Rix 2000) using an approximate, one dimensional treatment. We assume that the orbit 
remains circular, and solve for the coupled evolution of the disk surface density S(r, t) and 
secondary semi-major axis a. The governing equation for the surface density evolution is 
(Lin & Papaloizou 1986; Trilling et al. 1998) 



<9£ _ 1 d 
dt r dr 



(5) 



dr v ' (GM 1 ) 1 /2 

The first term on the right-hand side describes the usual diffusive evolution of an accretion 
disk. The second term describes how the disk responds to the torque from the secondary, 
which is approximated as a fixed function of radius A(r, a). We adopt 

fq 2 GMi ( r V 

A = ^— r < a 

2r \A p J 

fq 2 GM 1 ( a V 

= ^{aJ r>a (6) 

where / is a dimensionless normalization factor and A p is given by A p = max(/i, \r — a\). 
The rate of migration of the secondary is 

da ( a ^ 1/ V 47^ ^ r° ut • 

where the integral is taken over the entire disk. 

To fix /, and to help visualize the interaction, we have run high resolution (400 2 mesh 
points) simulations of the secondary-disk interaction using the ZEUS hydrodynamics code 
(Stone & Norman 1992). Figure 1 shows the results from a locally isothermal (c s = c s (r) 
only) calculation, which was evolved for 400 orbits of the secondary until a near steady-state 
was obtained. For both this run, with (h/r) ~ 0.07, and for simulations of cooler disks with 
(h/r) = 0.035, 0.02, the edge of the gap lies close to the 2:1 orbital resonance. A low level of 
accretion onto the secondary from the outer disk is observed. In our ID disk code described 
below, gaps of approximately the correct width are obtained using / = 10 -2 . 
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We complete the model by adopting a viscosity v — 1.5 x 10 20 (r/0.1 pc) 3//2 cm 2 s _1 , 
where the prefactor is chosen to match the viscosity of a disk with a = 10~ 2 and {h/r) = 0.01 
at 0.1 pc. The scaling with radius is model dependent and poorly known. We adopt a power- 
law with an index of 3/2 primarily for numerical convenience. 

3.1. Results for prompt mergers 

Figure 2 shows the evolution of a black hole binary with masses Mi = 5 x 10 s M Q and 
M2 = 10 7 M Q , calculated by numerically integrating equation (5) using an explicit method 
on a non-uniform grid. 800 mesh points were used to resolve the disk from r ms = 1.4 x 10~ 4 pc 
out to 3 pc. The disk has a zero-torque boundary condition at the inner edge, and an initial 
surface density profile corresponding to a constant accretion rate of 1 MQyr™ 1 out to 1 pc. 
Starting from a separation of 0.1 pc, a combination of disk-driven migration and inspiral due 
to gravitational radiation leads to merger within 2 x 10 7 yr. 

In Figure 3 we show the disk surface density profile during the migration and inspiral. 
There are two phases. In the first, the secondary migrates inwards within a gap under 
the action of disk torques, while the inner disk is partially depleted by accretion onto the 
primary black hole. The accretion rate onto the primary black hole in this first phase is 
somewhat reduced compared to the value it would have in the absence of a binary. This 
behavior can be easily understood - adding a binary companion to the disk merely increases 
the amount of angular momentum that must be transported outwards before mass can be 
accreted. Flux-limited samples of AGN may therefore select against finding close binaries. 
Disk-driven migration continues until angular momentum loss due to gravitational radiation 
first becomes important, at a separation of ~ 10~ 2 pc. Thereafter, the inspiralling black 
hole begins to sweep up the gas in the inner disk, forming a narrow inner spike in the disk 
surface density that is pushed inwards. The outer disk is unable to evolve on such a rapid 
timescale, so a large gap develops between the outer disk and the secondary. 

The final stages of the merger drive rapid accretion of the inner disk. Around 3 x 
10 5 M & of gas remain in the inner disk when the secondary reaches 10~ 3 pc, around 200 yr 
prior to merger. This implies an enormous accretion rate (formally exceeding 10 6 M Q yr _1 ) 
immediately prior to merger, and signals the breakdown of the thin-disk approximation. 
Instead of being radiated, the dissipated energy will instead go into thermal energy, puffing 
up the inner disk and forming a hot accretion flow (Begelman & Meier 1982). There are 
large uncertainties in the further evolution, but we speculate that the final merger is likely to 
occur within a common envelope formed from the disrupted inner disk. As in the common 
envelope evolution of binary stars (Taam, Bodenheimer & Ostriker 1978), strong outflows 
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are extremely likely, with characteristic velocities comparable to the orbital velocity of ~ 
10 4 kms" 1 . Strong outflows from galactic nuclei may therefore act as signposts to systems 
where binaries are either currently inspiralling, or have recently merged. 

3.2. Quiescent galactic nuclei from black hole binaries 

The timescale on which black holes form a hard binary, following a galactic merger, 
depends upon the stellar dynamics of the interaction between the binary and the stellar 
population. Gas could well flow into the merged nucleus on a shorter timescale (Hernquist 
& Mihos 1995). As a result, even if the merging galaxies contain substantial quantities of gas, 
it is not obvious that the resulting black hole binary will find itself in a gas-rich environment. 

Figure 2 also shows how the secondary migrates in disks that have the same viscosity 
as that considered previously, but initially lower disk masses and accretion rates. We ran 
one model with an initial accretion rate of 0.1 M yr _1 and a disk mass of 6 x 10 6 M Q , 
and a second with M = 1CT 2 M Q yr _1 and Mdi S k = 6 x 10 5 M & . Migration rates in these 
lower mass disks are substantially reduced, so it is possible that there could be a long-lived 
population of binaries embedded in low mass disks. The reduced level of accretion onto the 
binary, due to the gravitational torques, could then contribute to the surprisingly low levels 
of accretion inferred in most local galactic nuclei. This suppression of accretion, however, is 
only effective for cold gas. Hot gas, with low specific angular momentum, can surmount the 
angular momentum barrier and replenish small disks around the binary components (Bate & 
Bonncll 1997). Explaining why this accretion does not lead to large luminosities (Fabian & 
Rees 1995; Loewenstein et al. 2001) almost certainly involves the radiation physics (Narayan 
& Yi 1995), and possibly hydrodynamics (Blandford & Begelman 1999), of low M flows. 

4. Summary 

The ubiquity of supermassive black holes in the local Universe (Magorrian et al. 1998) 
suggests that black holes may also be common in the high redshift progenitors that merged 
to make up present day galaxies. Although a large black hole fraction at high-z is not proven 
(Menou, Haiman & Narayanan 2001), a high frequency of mergers between galaxies harboring 
black holes would have interesting implications for the accretion history of AGN (Haehnelt, 
Natarajan & Rees 1998; Haehnelt & Kauffmann 2000; Cavaliere & Vittorini 2000), and for 
the structure of galaxies themselves (Ravindranath, Ho & Filippenko 2002). 

In this letter, we examined the coupled evolution of a supermassive black hole binary 
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embedded in a gaseous accretion disk. First, we showed that, in a gas-rich environment, a 
binary with a separation of 0.1 pc merges within around 10 7 yr. Close binary quasars are 
therefore expected to be uncommon. Long lived black hole binaries are more likely to be 
found in apparently quiescent galactic nuclei at low redshift, where the disk accretion rate 
is lower, and where migration may be further slowed by inefficient disk angular momentum 
transport (Menou & Quataert 2001). Second, we examined the possible observational signa- 
tures of the merger. We argued that large enhancements to the accretion rate were unlikely 
during the phase of disk-driven migration. Forced accretion of any inner disk during the 
final, gravitational wave driven inspiral, however, is likely to lead to much of the inner disk 
being expelled in a high velocity quasi-spherical outflow. 

PJA thanks Yale for hospitality during the course of this work. Numerical simulations 
made use of the UK Astrophysical Fluids Facility. 
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Fig. 1. — The interaction between a low mass secondary black hole and the outer accretion 
disk, simulated at high resolution in two dimensions. The tightly wound spiral waves that 
mediate angular momentum exchange between the secondary and the disk are clearly seen. 
This run was for a mass ratio q = 0.01, disk thickness (h/r) ~ 0.07, and a constant kinematic 
viscosity v appropriate to a disk with a ~ 10~ 2 at the radius of the satellite. In general, 
similar waves would also propagate through an inner accretion disk, which was not included 
in this calculation. 
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Fig. 2. — Migration of the secondary black hole in disks with initial accretion rates of 
1 M yr _1 (solid line), 0.1 M yr _1 (short-dashed line) and 10~ 2 M yr -1 (long-dashed line). 
Binaries embedded within relatively massive disks are quickly driven into the gravitational 
wave dominated regime and merge promptly. 
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Fig. 3. — Disk surface density during black hole migration. The black hole migrates inwards 
within a gap during the disk dominated phase (upper panel, curves show S at three epochs). 
Once orbital decay due to gravitational radiation dominates, the inner disk is swept up by 
the inspiralling hole and forms a surface density spike (lower panel). The outer disk, which 
is unable to evolve on such a short timescale, maintains a fixed surface density profile. 



